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We extend recent ab initio calculations of the electronic band structure and the phonon dispersion 
relations of rhombohedral GeTe to calculations of the density of phonon states and the temperature 
dependent specific heat. The results are compared with measurements of the specific heat. It is 
discovered that the specific heat depends on hole concentration, not only in the very low temperature 
region (Sommerfeld term) but also at the maximum of Cp/T 3 (around 16 K). To explain this 
phenomenon, we have performed ab initio lattice dynamical calculations for GeTe rendered metallic 
through the presence of a heavy hole concentration (p ~ 2x 10 21 cm -3 ). They account for the 
increase observed in the maximum of Cp/T 3 . 



PACS numbers: 63.20.D-, 65.40.Ba 

I. INTRODUCTION 

GeTe is a material with 10 valence electrons per prim- 
itive unit cell (PC) and is thus related to the lead chalco- 
genides PbX (X = S, Se, Te) and to the semimetals As, 
Sb, and Bi. Like PbX, GeTe is a semiconductor but crys- 
tallizes in the rock salt structure only at temperatures 
above - 700Ki Below - 650K, GeTe distorts through 
an extension of the cube diagonal into a rhombohedral 
structure^ similar to that of the semimetals As, Sb and 
Bi, except that the two atoms per PC are not equal and 
hence its bonding is polar. Compatible with the rhom- 
bohedral symmetry represented by a strain along [111], 
there are two additional independent distortion parame- 
ters, the rhombohedral angle a, whose distortion from 
the 60° of the cubic phase brings it to «58°, and a 
change in the Ge-Te sublattice separation aor\/3 with 
t - 0.03J- 3 . 

This distortion converts GeTe into a ferroelectric with 
T c ~ 650 K, given the polar nature of the bond. An- 
other interesting property of GeTe is its native p-type 
doping, with hole concentrations typically higher than 
5x 10 19 cm~ 3 . This fact is attributed to the presence of 
Ge vacancies 4^ 

The metallic nature of GeTe with high hole concentra- 
tion (p > 8.3xl0 20 ) leads to superconductivity at tem- 
peratures T c below ~ 0.3 K£ Finegold^ examined the 
heat capacity of a GeTe sample with p ~ 1.1 xlO 21 cm -3 
and T c ~ 0.25 K in the region around T c and observed the 
peak related to superconducting behavior: application of 
a magnetic field of 500 Oe completely wiped out super- 
conductivity and the related peak, thus allowing him to 
determine, for his sample, a Sommerfeld term 7 = 1.32 
mJ/moleK 2 . The specific heat measurements reported 
here were performed from 4 to 250 K, and therefore do 
not address the superconducting behavior of GeTe. They 
allow us, however, to determine 7 and the maximum in 
Cp/T 3 ( at ~17 K) which signals deviations from Debye's 



T 3 law. We find this maximum to be 5.5% higher than 
that obtained from the data of Bevolo et alA, a fact that 
suggests a decrease in phonon frequencies induced by the 
presence of hole doping whose concentration is, in our 
sample, ~25 times higher than in that of Bevolo et al.. 

In order to test this hypothesis we have performed ab 
initio calculations of the lattice dynamics of a-GeTe simi- 
lar to those reported earlier for PbS^ and for PbS, PbSe, 
and PbTe 9 -. Comparison of theoretical -obtained with- 
out spin-orbit (SO) coupling- and experimental values of 
Cp/T 3 presented in Ref. Q suggested that there may be 
a significant contribution (~ 20%) of the SO interaction 
to the value of C p /T 3 at its maximum (note that in the 
region of the maximum C p ~ C v &) 

We therefore repeated in Ref. d lattice dynamical cal- 
culations for PbS, PbSe and PbTe, this time with and 
without SO interaction. We were able to prove that the 
lack of SO interaction, especially at the Pb ions, was 
largely responsible for the discrepancy just mentioned. 
In the calculations reported here for cv-GeTe we have 
checked the effect of SO interaction and, not unexpect- 
edly (heavy Pb ions are not present), we found it to 
be rather small (~1%); it can thus be neglected in the 
remaining of this work. We did find, however, for the 
calculated values of C p /T 3 an increase in its maximum 
between an undoped sample and one with p ~ 2.1xl0 21 
cm -3 of ~6%, which largely explains the difference found 
between our sample and that of Bevolo et al. (3.5%). Our 
calculations of the lattice dispersion relations for p ~ 
2.1xl0 21 cm -3 and for undoped samples show indeed 
that the phonon frequencies decrease with doping. This 
effect is very large around the T-point of the Brillouin 
zone (BZ) due to the screening of the ionic charge by the 
holes. However, it amounts to a few percent throughout 
most of the BZ. 

As a by-product we have also calculated the phonon 
density of states (DOS) and its projection on each of the 
constituent ions. We have also calculated the density of 
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two-phonon states which applies to optical spectroscopics 
involving two phonons: second order Raman scattering 
and infrared absorption. It is hoped that the availabil- 
ity of these DOS calculations will encourage such optical 
measurements in this interesting material . 

II. THEORETICAL DETAILS 

All the calculations have been performed using plane 
waves and norm-conserving pseudopotcntials, as imple- 
mented in the ABINIT code i n i 12 The dynamical prop- 
erties have been evaluated within the density-functional 
perturbation theory 13 i 14 i 15 We employed Hartwigsen- 
Goedecker-Hutter (HGH) pseudopotentials^ generated 
including spin-orbit coupling, within the local den- 
sity approximation (LDA) adopting the Teter-Pade 
parameterization^ Even though HGH pseudopotcntials 
are known as being relatively hard, the properties inves- 
tigated in this work are well converged when including 
a plane wave basis up to a kinetic energy cutoff equal 
to 15 Ha. The Brillouin zone integration was performed 
using special k-points sampled within the Monkhorst- 
Pack schemed We found that a mesh of 12 x 12 x 12 
k-points was required to describe well the structural and 
vibrational properties. Further details on the theoretical 
calculations arc mostly given in Ref. except those for 
the heavily p-doped material which were only performed 
without SO interaction. 

To mimic the effect of Gc vacancies on the calculated 
force constants, and consequently on the heat capacity, 
the number of total electrons was reduced by a certain 
amount Se. The value of Se was chosen to be equivalent to 
the concentration of the free holes p. In all calculations of 
hole doped cases, a background of homogeneous negative 
charge was imposed to achieve charge neutrality 

Before calculating the dynamical matrix the structural 
parameters clq , a, and r were optimized. In Table U 
we show the calculated structural parameters with and 
without SO coupling. The effect of the SO interaction 
on the structural parameters is found to be negligible. 



TABLE I: Calculated structural parameters of GeTe. The lat- 
tice parameter ao (in A ), the angle a (in deg), the deviation 
from of the Ge sublattice from the 0.5 sublattice position r, 
and the volume f2 = (ajj/4) sin 2 a (in A ). 





ao 


a 


T 


n 


with SO 


5.897 


88.96 


0.0237 


51.26 


without SO 


5.894 


88.96 


0.0237 


51.19 


Theory^ 


5.886 


89.24 


0.0217 


50.96 


Exp -2 


5.996 


88.18 


0.026 


53.84 


Exp -22 


5.98 


88.35 


0.0248 


53.31 



Our results are globally in good agreement with previ- 
ously reported ab initio results.—. Moreover the devia- 
tion between the calculated structural parameters(ao,a) 
and experiment is less than 2%. We performed the calcu- 
lations using the parameters obtained without spin-orbit 
coupling for all cases including those of doped samples. 



III. CALCULATED DISPERSION RELATIONS 
AND PHONON DENSITY OF STATES 

In Fig. [T]we show the phonon dispersion relations cal- 
culated for undoped GeTe (p = 0) and for GeTe contain- 
ing p = 2.1xl0 21 holes/cm 3 . The undoped material, 
to the best of our knowledge, has never been prepared. 
We calculate it as a reference point to assess the effects 
of p-type doping on the dispersion relations. For p = 
we have performed calculations without and with SO in- 
teraction. This interaction has a very weak effect, only 
sizeable near the T-point of the BZ: the frequency of the 
LO phonons near T is lowered by 1.3 cm -1 when the 
SO interaction is included (see Table llT|) . This effect can 
be compared with that found for PbTe^ for the equiva- 
lent phonons, which lies around 5 cm -1 . Assuming that 
this effect is proportional to the sum of SO splittings of 
the valence electrons of cation and anion^ (for values of 
atomic SO splittings see Ref l2ll ) we predict, from the 
value of 5 cm -1 calculated for PbTe, 1.4 cm -1 for GeTe, 
compatible with the value of 1.3 cm -1 obtained in the 
calculations. The softenings induced by SO coupling for 
other phonons in the BZ are even smaller and thus will 
be henceforth neglected. 




FIG. 1: (color online) Calculated phonon dispersion relations 
of GeTe along selected symmetry directions in the BZ. The 
(black) solid line denotes the results for stoichiometric un- 
doped GeTe, the (red) solid line denotes the results including 
spin-orbit coupling and the (blue) solid line was obtained for 
doped GeTe with p = 2xl0 21 holes/cm 3 without spin-orbit 
coupling. 

In Fig. Q]we can also appreciate the effect of hole dop- 
ing on the dispersion relations. Without doping, the TO 
phonons are split by as much as 45 cm -1 on account 
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of the electrostatic field induced by the large transverse 
ionic charges^. This effect is a non-analytic function 
of wave vector k : for fc « the LO and TO frequen- 
cies depend on the direction of k. The corresponding 
discontinuities disappear in the doped material, as the 
transverse charges are screened by the free holes: like in 
the cases of As, Sb and B i 22 ' 23 the optical phonons of the 
doped GeTe are split by the rhombohcdral field into a sin- 
glet (LO) and a doublet (TO), the splitting amounting 
to 42 cm -1 , in agreement with measurements by Ra- 
man spectroscopy^ and with our experimental results 
(see Table [TTI and section IV- A). The individual calcu- 
lated frequencies are, however, smaller than experimen- 
tal results. The differences between the absolute values 
of calculated and measured frequencies is probably due 
to the employed local density and pscudopotential ap- 
proximations, or, to temperature effects resulting from 
anharmonic terms which we do not include in the present 
calculations. Similar behavior has been noticed for the 
lead chalcogenides where the calculated LO frequency of 
PbTe is 28 cm -1 lower than the measured one, see Fig. 
3 of Ref. ©. 

We should mention at this point that similar soften- 
ings of the phonons of semiconductors upon heavy dop- 
ing have been investigated earlier, both experimentally 
and theoretically. Most recently the case of boron doped 
diamond has been considered, in connection with the su- 
perconductivity observed in this materiaL 25 i 26 For a hole 
concentration similar to the one used here Boeri et al. 
find for diamond a softening of the optical frequency at 
T of ~20%. From Table [TTI we estimate a doping induced 
softening of 15% for the LO phonon and 30% for the TO 
phonons at T. It has been pointed out that the strong 
phonon softening, in the case of diamond, is related to 
a strong hole-phonon interaction, probably also respon- 
sible for superconductivity^ Both, ferroelectricity and 
superconductivity, point to a strong carricr-phonon in- 
teraction in GeTe. Phonon softenings upon doping with 
boron have also been observed in Sij 28 i 29 

We have calculated the phonon DOS and its projec- 
tions on the Ge and Te atoms using the same procedure 
as m Refs. @and!. 

The one-phonon DOS is shown in Figj2ja) for three 



TABLE II: Phonon frequencies at the zone center (in cm x ) 
calculated for a q vector || to trigonal axis. 
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p = (no spin-orbit) 


73 


153 


p — (with spin-orbit) 


72 


152 


p = 2.1 x 10 21 holes/cm 3 


73 


115 


Exp (present work) 


88 


123 


Exp2i 


80 


122 
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FIG. 2: (color online) (a) One-phonon DOS of a-GeTe calcu- 
lated for an intrinsic sample (solid lines) without SO interac- 
tion, (dotted line) with SO interaction included, and (dashed 
line) for a sample with p = 2.1xl0 21 holes/cm 3 , (b) Pro- 
jection of the DOS curves of Fig. 2a on Ge (red) and on Te 
(green) atoms. 



different cases: undoped GeTe without SO interaction 
and with SO interaction included (note that the differ- 
ence is hardly visible and, as already mentioned above, 
can be neglected) and GeTe with p = 2.1x 10 21 holes/cm 3 
(dashed line, no SO interaction). The presence of holes at 
the top of the valence band lowers all of the three bands 
in Fig. [D^a), as corresponds to the softening of the dis- 
persion relations in Fig. [TJ Figure^b) displays the three 
DOS of Fig. HJa) projected on the two constituent atoms 
of a-GeTe. In Figure [2] three bands are observed. They 
correspond to the acoustic phonons (from to 110 cm" 1 ), 
the TO, and the LO phonons, respectively. Whereas in 
PbX (X=S, Se, Te)£ the TA band is almost exclusively 
Pb-likc, Fig. [2^b) shows that for GeTe this band con- 
tains a strong mixture of Te and Ge vibrations (~60- 
40%) whereby the heavier ion dominates, as expected for 
acoustic modes. The intermediate (TO) band contains a 
nearly equal mixture of vibrations of both atoms, with a 
slight predominance of the lighter one (~ 90-10 %). The 
LO-band is clearly dominated by Ge vibrations (~30- 
70%), as expected. All bands are lowered in the heavily 
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p-dopcd material with respect to the undoped ones, as 
expected from the dispersion relations in Fig. [TJ This 
lowering results in an increase in the maximum of C„. p 
vs. T found at ^16K, as will be discussed in Sect. HVl 



IV. EXPERIMENTAL PROCEDURE AND 
RESULTS 

A. Samples and Sample Characterization 

Most of the measurements reported here were per- 
formed on wafers of a coarse polycrystalline ingot piece 
grown by the Bridgman technique in a quartz ampoule^ 
The ingot diameter was about 10 mm. We cut from it 
several wafers perpendicular to the growth axis, about 
0.8 mm in thickness. The wafer measured by us will be 
labeled GeTe-CG. We also measured the heat capacity of 
commercial samples (99.999% purity, Alfa Aesar, Karl- 
sruhe, Germany) consisting of small crystalline chips, 
about 1 mm in size. Several chips weighing a total of 
~23 mg were measured (they are labeled as Alfa-GeTe 
below). The results are compared with the data pub- 
lished by Bevolo et alA. 

The wafers were characterized by Hall effect measure- 
ments using the Van der Pauw technique^. No such 
measurements were possible for the granular Alfa Aesar 
sample. In this case, the hole concentration p was es- 
timated from the linear-in-T term of the specific heat 
(Sommerfeld term, equal to jT) by comparing it with 
that reported in Ref. H for p = 8xl0 19 holes/cm 3 . For 
this purpose we used the relation^ 2 .: 

7 cx m*p 1/3 (1) 

where m* is the effective band mass of the holes. 
We also estimated p in the CG2-GeTe sample from the 
plasma minimum in the ir reflection spectrum by com- 
parison with similar data in Ref. HU. We found for this 
sample p= 1.5x 10 21 cm" 3 . The values of p and 7 for the 
samples under consideration and the estimated values of 
p arc listed in Table UTTl 

Table IIHI displays the hole concentration of our two 
samples and that given by Bevolo^ as determined by Hall 
measurements (ph) together with the concentration es- 
timated from the ir reflectivity (7 op ) and those obtained 
from the measured values of 7. The latter were estimated 
by scaling the values of Bevolo et al. using Eq. (TTJ . In 
doing so, we have taken into account the fact that the 
effective mass m* increases with increasing p due to non- 
parabolicity of the bands. Since p 7 is proportional to 
(m*) , its estimated value depends very critically on this 
nonparabolic increase, which is not accurately known at 
present. We have taken for this increase the values re- 
ported in Ref. [34|: 24% for the Alfa Aesar sample and 
32% for the CG sample with respect to Bcvolo's result. 
In spite of some variations, the values of p listed in this 



table are adequate for the subsequent discussion of our 
heat capacity measurements. 

For the sake of completeness we also measured the Ra- 
man spectrum of the LO and TO phonons in our CG- 
GeTe sample. We found at 293 K wlo = 123 cm" 1 and 
wto = 88 cm" 1 , in good agreement with the values re- 
ported in Ref. |35j for the same temperature. 

B. Temperature dependence of the heat capacity 

We display in Fig. [3] the heat capacity of GcTe mea- 
sured below 5.5 K using the standard Sommerfeld plot 
{C p /T vs. T 2 ) : The intercept of a linear fit with the 
ordinate yields the value of 7 as listed in Table UTTl 




T 2 (K 2 ) 



FIG. 3: (color online) Sommerfeld plots of the heat capacities 
of our CG and Alfa samples of a-GeTe together with that for 
Bevolo's sampled The 7 values obtained from the intercept 
with the vertical axis are listed in Table [TTT1 

Figure 2] shows the temperature dependence of C p /T 3 
measured for our two samples in the 2.5 - 30K temper- 
ature range, together with the results of Bevolo et al.A 
The increase in C p /T 3 seen below 5 K corresponds to the 
Sommerfeld term shown in Fig. [3] The most remarkable 
feature of Fig. 2] is the monotonic increase in Cp/T 3 ob- 
served with increasing p, which can be assigned to the 
frequency softening (cf. Fig. [TJ and the corresponding 
down-shift of the DOS (cf. Fig. HJa)). Notice also that 
our two samples exhibit this maximum at 15.8 K whereas 
Bevolo's maximum is found at 16.6 K. This shift can be 
assigned to the fact that the latter sample has a hole 
concentration over an order of magnitude lower than the 
former. 

A simple relationship has been proposed earlier (Ref. 
0H) to relate the temperature of the maximum in C p /T 3 
to the phonon DOS. One first divides the temperature 
by 1.43 so as to convert T (in K) into cm" 1 . Then one 
divides the obtained frequencies by 6.2. This procedure 
leads to the frequencies 68.5 cm" 1 for our samples and 
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TABLE III: Values of the hole concentration of our samples in cm 3 as determined by various methods (see text). Also, values 
of the Sommerfeld coefficients 7 in mJ/molK 2 and f3 = 2xl2/5-7r 4 R/#D e bye in mJ/molK 4 are given as well as the value of 
Cp/T 3 at the maximum in mJ/molK 4 . A blank signifies that a value is not available. The values from Ref. are also cited for 
comparison. 
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Cp/2~niax 


Pi 


7o P 


Bevolo 


8xl0 19 


0.554(3) 


0.307(1) 


0.640(2) 


8xl0 19 




CGI 


1.6xl0 21 


1.66(9) 


0.380(9) 


0.663(2) 


1.2xl0 21 




CG2 




1.57(2) 


0.368(3) 


0.663(2) 


1 xlO 21 


1.5 xlO 21 


CG3 




1.19(2) 


0.339(2) 


0.658(2) 


4 xlO 20 




Alfa 




1.45(4) 


0.358(3) 


0.656(2) 


8 xlO 20 





700 I , , , 1 , 1 , 1 r 




FIG. 4: (color online) Measured temperature dependence of 
Cp/T 3 for the three samples of a-GeTe under consideration. 
Details for other samples are given in Table ITTT1 Note the in- 
crease of Cp/T 3 with increasing hole concentration (see Table 

EE). 



FIG. 5: (color online) Calculated temperature dependence of 
C v /T 3 for an undoped a-GeTe sample (with and without SO 
splitting), and for a sample with p = 0.5 x 10 21 , 1.0 x 10 21 , 
and 2.1 x 10 21 cm~ 3 (without SO splitting). 



72 cm -1 for Bevolo's sample.— The corresponding shift 
of 3.5 cm -1 agrees with that observed for the first TA 
peak in the DOS. The position of this peak, 70 cm -1 
also agrees with the average position estimated by the 
procedure followed above. 

Figure [5] displays the temperature dependence of 
C v /T 3 (~ Cp/T 3 ) calculated for several different concen- 
trations (p=0, 0.5xl0 21 , l.OxlO 21 , and 2.1xl0 21 cm- 3 ). 
Notice that the effect of SO coupling, as implemented 
for the undoped sample, is small and we have not imple- 
mented it for doped GeTc. Hence, in the latter case we 
should compare it with the undoped calculation without 
SO interaction. We did not notice much of a change for 
the light doping case (p=0, 0.5xl0 21 cm -3 ). However, 
as expected, the increase in doping concentration leads 
to a monotonic increase in the maximum of C v /T 3 . As 
seen in Figure [3 we found an increase in the maximum of 
0.2 and 0.6% for p^l.OxlO 21 and 2.1 x 10 21 holes/cm 3 , 
respectively. The results agree rather well with the 0.5% 
measured for the CG sample, especially when one con- 



siders that this sample has a doping of about 1.5 xlO 21 
cm~ 3 . The shift in the calculated position of the maxi- 
mum with doping is also close to the measured one. 

It is of some interest to compare the small effect of the 
SO interaction shown in Fig. [5] with that calculated for 
PbTe. In Ref. we proposed a perturbation expansion 
of the SO effect as a function of the strength of the SO 
splittings of the p valence electrons of Pb and Te. We 
use the same expression replacing the SO splitting of Pb 
(L26 eV, cf. Ref. M) by that of Ge (0.22 eV, cf. Ref. 
l2lh and predict an effect of the SO coupling of 4% for 
the maximum of C v /T 3 in a-GeTe (as opposed to 20% 
for PbTe). The effect observed in Fig. lis about 2%, 
somewhat smaller than the rough estimate carried over 
from the PbTe. 

We conclude by displaying in Fig. [5] the temperature 
dependence of C p as measured by us, Bevolo et alA and 
Zhdano'v 3 ^, compared with our ab initio calculations of 
C v up to 300 K. Our measurements up to 150 K, from 
which the Sommerfeld term has been subtracted, agree 
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where a v (T) is the temperature dependent coefficient 
of the volume thermal expansion, B is the bulk modulus 
and V mo \ is the molar volume. 

We have evaluated Eq.jH ) at 300 K using 
a.=5.60xl0- 5 K- 1 (Ref. M), B = 49.9 GPa (Refs. 
T0II20] ) and a molar volume of 32.3 cm 3 . We find a 
thermal expansion contribution to C p - C v of ~ 1.5 
J/molK at 300 K, which is by a factor of ^2.5 smaller 
than the value required to bring Zhdanov's data (C p ) 
to agree with the calculation (C v ). In order to clarify 
this matter, measurements should be repeated in the 
temperature range 150 - 300 K. 



FIG. 6: (color online) Heat capacity divided by T 3 measured 
for samples (Bevolo and our data CGI) compared with calcu- 
lations withp=0 and p=2.1x 10 21 cm -3 without SO coupling. 




300 



7(K) 



FIG. 7: (color online) Measured heat capacity of several a- 
GeTe samples in the 3-300 K range. The black solid line 
represents our ab initio calculations. The linear dependence 
in Zhdanov's data for T < 150 K corresponds to the thermal 
expansion effect. 

within error, with those of Zhdanov. Between 150 and 
300 K, Zhdanov's data lie above the calculations. In 
order to ascertain whether this discrepancy is due to the 
difference between C p and C v we use the expression: 

C p (T) - C v (T) =a 2 v (T)-B- V moi ■ T, (2) 



V. CONCLUSIONS 

We have presented a detailed experimental and theo- 
retical investigation of lattice dynamics and heat capac- 
ity of a-GeTe. In contrast to previously reported results 
for other IV- VI materials^, we found that the inclusion 
of spin orbit coupling has a small effect on the calcu- 
lated phonon frequencies and consequently on the heat 
capacity. On the other hand, we found that ignoring 
non-stoichiometry effects of GeTe, leads to discrepancies 
between the calculated and the measured C p /T 3 . These 
can be resolved by taking into account the existence of 
free holes which lead to an increase of C p /T 3 . Such in- 
crease was found to be due to phonon frequency soften- 
ing and the corresponding down shift of the density of 
phonon states. 
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